AJ, in press (Aug, 2004) 



A Survey of z > 5.7 Quasars in the Sloan Digital Sky Survey III: Discovery of 

Five Additional Quasars 1 

Xiaohui Fan 2,3 , Joseph F. Hennawi 3,4 ', Gordon T. Richards 4 , Michael A. Strauss 4 , Donald P. 
Schneider 5 , Jennifer L. Donley 2 , Jason E. Young 2 , James Annis 6 , Huan Lin 6 , Hubert Lampeitl 6 , 
Robert H. Lupton 4 , James E. Gunn 4 , Gillan R. Knapp 4 , W. N. Brandt 5 , Scott Anderson 6 , Neta 
A. Bahcall 4 , Jon Brinkmann 7 , Robert J. Brunner 8 , Masataka Fukugita 7 , Alexander S. Szalay 9 , 

Gyula P. Szokoly 10 , Donald G. York 11 

ABSTRACT 

We present the discovery of five new quasars at z > 5.7, selected from the multicolor 
imaging data of the Sloan Digital Sky Survey (SDSS). Three of them, at redshifts 5.93, 



1 Based on observations obtained with the Sloan Digital Sky Survey, and with the Apache Point Observatory 
3.5-meter telescope, which is owned and operated by the Astrophysical Research Consortium; and with the MMT 
Observatory, a joint facility of the University of Arizona and the Smithsonian Institution, with the Univeristy of 
Arizona 2.3-meter Bok Telescope, with the Kitt Peak National Observatory 4-meter Mayall Telescope, with the 
6.5-meter Landon Clay Telescope at the Las Campanas Observatory, a collaboration between the Observatories 
of the Carnegie Institution of Washington, University of Arizona, Harvard University, University of Michigan, 
and Massachusetts Institute of Technology , and with the Hobby-Eberly Telescope, which is a joint project of 
the University of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximillians- 
Universitat Miinchen, and Georg-August-Universitat Gottingen. 

2 Steward Observatory, The University of Arizona, Tucson, AZ 85721 

3 Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement 
with the National Science Foundation. 

4 Princeton University Observatory, Princeton, NJ 08544 

5 Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802 
6 Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510 
6 University of Washington, Department of Astronomy, Box 351580, Seattle, WA 98195 
7 Apache Point Observatory, P. O. Box 59, Sunspot, NM 88349-0059 

institute for Cosmic Ray Research, University of Tokyo, Midori, Tanashi, Tokyo 188-8502, Japan 
8 Dept. of Astronomy & NCSA, University of Illinois, 1002 W. Green Street, Urbana, IL 61801 
9 Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218, USA 
10 Max-Planck-Institut fur extraterrestrische Physik, Postfach 1312, 85741 Garching, Germany 
11 University of Chicago, Astronomy & Astrophysics Center, 5640 S. Ellis Ave., Chicago, IL 60637 



- 2 - 



6.07, and 6.22, were selected from ~ 1700 deg 2 of new SDSS Main Survey imaging 
in the Northern Galactic Cap. An additional quasar, at redshift 5.85, was discovered 
by coadding the data obtained in the Fall Equatorial Stripe in the SDSS Southern 
Survey Region. The fifth object, at redshift 5.80, is selected from a non-standard 
SDSS scan in the Southern Galactic Cap outside the Main Survey area. The spectrum 
of SDSS J162331.81+311200.5 (z = 6.22) shows a complete Gunn-Peterson trough at 
z a bs > 5.95, similar to the troughs detected in other three z ^ 6.2 quasars known. 
We present a composite spectrum of the z > 5.7 quasars discovered in the SDSS to 
date. The average emission line and continuum properties of z ~ 6 quasars exhibit no 
significant evolution compared to those at low redshift. Using a complete sample of 
nine z > 5.7 quasars, we find that the density of quasars with M1450 < — 26.7 at z ~ 6 
is (6 ± 2) x 10~ 10 Mpc~ 3 (H = 65 km s^Mpc -1 , Q = 0.35 and A = 0.65), consistent 
with our previous estimates. The luminosity distribution of the sample is fit with a 
power law luminosity function ^f(L) oc x- 3 - 2 ^ - 7 ^ somewhat steeper than but consistent 
with our previous estimates. 

Subject headings: quasars: general; quasars: emission line; quasars: absorption lines 

1. Introduction 

This paper is the third in a series presenting i-dropout z ^ 5.7 quasars selected from the 
multicolor imaging data of the Sloan Digital Sky Survey (SDSS; York et al. 2000, Stoughton et 
al. 2002, Abazajian et al. 2003, 2004). In Fan et al. (2000) and in the first two papers of this 
series (Fan et al. 2001c, Paper I, Fan et al. 2003, Paper II), we presented the discovery of seven 
luminous quasars at z = 5.74 — 6.42, selected from ~ 2900 deg 2 of SDSS imaging in the Northern 
Galactic Cap. In this paper, we describe the discovery of five quasars at z = 5.80, 5.85, 5.93, 6.07 
and 6.22, respectively. The scientific objectives, photometric data reduction, target selection and 
follow-up observation procedures are described in detail in Paper I. Three of the new quasars 
were selected using the same color selection procedures described in Paper II, as outlined briefly 
in §2.1. One object was selected using the co-added catalog from multi-epoch imaging of the Fall 
Equatorial Stripe in the SDSS Southern Survey region. The co-addition procedures are described 
in §2.2. The final object was selected in a non-standard SDSS scan in the Southern Galactic 
Cap outside the main survey area. We present the spectroscopic follow-up observations of the 
i-dropout candidates and the photometric and spectroscopic properties of the newly-discovered 
quasars in §3. Combining the quasars in this paper and those presented in Fan et al. (2002) and 
in Papers I and II, we construct the composite spectrum of z ~ 6 quasars, and compare it with 
the average quasar spectrum at low redshift (§4). Finally, we update our estimate of the evolution 
of the high-redshift quasar luminosity function in §5. 

Following the previous papers in this series, we use two cosmologies to present our results: (1) 
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H = 50 km s- 1 Mpc- 1 , Q A = and Q M = 1 (fi-model); (2) H = 65 km s" 1 Mpc" 1 , ft A = 0.65 
and SIm = 0.35 (A-model). Coincidentally, the luminosity and number densities measured 
under our A-model are within 1% of those in the WMAP (Spergel et al. 2003) cosmology with 
H = 71 km s~ l Mpc' 1 , Q A = 0.73 and fl M = 0.27. 

2. Candidate Selection and Identification 

The Sloan Digital Sky Survey is using a dedicated 2.5m telescope and a large format CCD 
camera (Gunn et al. 1998) at the Apache Point Observatory in New Mexico to obtain images in 
five broad bands (u, g, r, i and z, centered at 3551, 4686, 6166, 7480 and 8932 A, respectively; 
Fukugita et al. 1996) of high Galactic latitude sky in the Northern Galactic Cap. About 6000 deg 2 
of sky have been imaged at the time of this writing (Mar 2004). The imaging data are processed 
by the astrometric pipeline (Pier et al. 2003) and photometric pipeline (Lupton et al. 2001), and 
are photometrically calibrated to a standard star network (Smith et al. 2002, see also Hogg et al. 
2001). In addition to the Main Survey, the SDSS also obtains multi-epoch imaging in ~ 270 deg 2 
along the Celestial Equator in the Southern Galactic Cap (the SDSS Southern Survey, York et al. 
2000). At the time of this writing, the Southern Survey region has been imaged between 5 and 15 
times, depending on the RA of the field. The multi-epoch data are used to study variable objects 
and are co-added to reach fainter limiting magnitudes. We will use imaging from both the Main 
Survey and the Southern Survey to search for quasar candidates. 

2.1. Selection in the Main Survey Area 

The SDSS quasar target selection pipeline (Richards et al. 2002) selects only quasar 
candidates at z <J 5.5. The highest-redshift object discovered in the SDSS spectroscopic survey is 
at z = 5.4 (Anderson et al. 2001). At higher redshift, quasars become i-dropout objects; their 
selections require additional follow-up observations. Papers I and II present the results from a 
survey of i-dropout candidates selected from ~ 2900 deg 2 of SDSS Main Survey imaging carried 
out in the Springs of 2000, 2001 and 2002. In Spring 2003, we searched for i-dropout quasar 
candidates in 37 new SDSS imaging runs. These imaging data were taken between 10 October 
2002 (Run 3358) and 29 April 2003 (Run 3919). We used the same criterion to decide which 
photometric runs to include in the i-dropout survey as in Papers I and II: the z band image 
quality, measured by the psfWidth parameter (= 1.06 FWHM for a Gaussian profile) in the fourth 
column of the SDSS camera, should be better than 1.8". The median seeing in the i and z bands 
is ~ 1.4" for the entire survey area (Abazajian et al. 2003). There is overlap between adjacent 
SDSS strips and stripes (York et al. 2000), meaning that these 37 new runs overlap somewhat 
with the area covered in Papers I and II. Taking these overlaps into account, we find that the 
total new area of the sky covered by these runs is 1708 deg 2 , bringing the combined sky coverage 
of Papers I, II and this paper to 4578 deg 2 . 



-4 - 



We applied the same color selection criteria as in Paper II (see Figure 1 and 2 in Paper II) 
to the new SDSS imaging data to selection z > 5.7 quasar candidates. A total of 80 i-dropout 
candidates were selected in the main survey area. The photometric and spectroscopic follow-up 
observations were carried out over a number of nights between December 2002 and June 2003. 
Independent z photometry was carried out using the Seaver Prototype Imaging camera (SPICAM) 
in the SDSS z filter on the ARC 3.5m telescope at the Apache Point Observatory. J photometry 
was carried out using the 256 x 256 NICMOS imager on Steward Observatory's 2.3m Bok Telescope 
at Kitt Peak, and using the GRIM II instrument (the near infrared GRIsm spectrometer and 
IMager), also on the ARC 3.5m. The spectroscopic follow-up observations were obtained using 
the Red-Channel Spectrograph on the MMT 6.5m telescope on Mt. Hopkins, the Double Imaging 
Spectrograph (DIS) on the ARC 3.5m, the Multi-Aperture Red Spectrograph (MARS) on the 
4-m telescope on Kitt Peak, and the Low Resolution Spectrograph (LRS, Hill et al. 1998) on the 
Hobby-Eberly Telescope. 

2.2. Selection in the SDSS Southern Survey Region 

The SDSS Main Survey imaging consists of single epoch observations with exposure times 
of 54.1 seconds. For the average seeing conditions (FWHM ~ 1-4"), the 5-cr limiting magnitudes 
in the i and z bands are 22.5 and 20.5, respectively (Stoughton et al. 2002, Abazajian et al. 
2003). At this depth, the SDSS Main Survey only allows selection of i-dropout candidates at 
z < 20.2, and therefore is only sensitive to the most luminous quasars at z ~ 6 (Mb <S —27). 
The multi-epoch imaging obtained in the SDSS Southern Survey area will eventually provide 
photometry ~ 1.5 magnitudes deeper than the main survey over ~ 270 deg 2 , enabling selection of 
much fainter quasar candidates. By the end of 2003, the SDSS Southern Survey Region had been 
scanned 5-15 times, resulting in a dataset that goes more than one magnitude deeper than the 
Main Survey. 

In Fall 2002, we used a preliminary co-added catalog of ~ 5 epoch imaging to select z ~ 6 
quasar candidates. The co- added catalog was generated by matching sources detected in more 
than one SDSS run within a radius of 1". The average flux is calculated from the SDSS asinh 
magnitudes (Lupton, Gunn Sz Szalay 1999), weighted by the inverse of the flux variances. Note 
that when co-adding at the catalog level, an object is required to be detected in the individual 
runs. Therefore, while the co-addition improves the S/N of faint sources, allowing us to select 
candidates all the way down to the 6-cr detection limit, it does not allow measurement of objects 
fainter than the detection limit of single SDSS exposures. Therefore, we can only select candidates 
at z < 20.5 — 20.7. We searched for candidates from ~ 100 deg 2 of co-added catalogs. The 
color selection criteria in the Southern Survey region are the same as those used in the Main 
survey area. A number of faint i-dropout candidates were observed using Magellan II Clay 
Telescope and the Boiler and Chivens (B&C) spectrograph on 16 Oct 2002. One quasar, SDSS 



- 5 - 



J000552.34-000655.8 14 (hereafter SDSS J0005-0006, z AB = 20.5) was discovered at a redshift of 
5.85. 



2.3. Selection in the Constant-Longitude Scans 

The SDSS has obtained some imaging data outside the survey boundary (York et al. 2000; 
Finkbeiner et al. 2004), including a few strips in the Southern Galactic Cap that extend to low 
Galactic latitude. We have also selected i-dropout candidates in these runs. One new quasar, 
SDSS J000239.39+255034.8 (z = 5.80), was selected from Run 4152, a constant-longitude scan (at 
/ = 110°), observed on 29 Sep 2003. The sample in these scans which go outside the main survey 
area is not yet complete, and so we do not include this object in our complete sample. 

3. Discovery of Five New Quasars at z > 5.7 

Among the 80 i-dropout candidates candidates in the Main Survey area, 15 are false z 
band only detections which are most likely cosmic rays; 55 are M or L dwarfs (mostly classified 
photometrically based on their red z* — J colors); and 7 are likely T dwarfs. Several objects 
still lack proper infrared spectroscopy, so the T dwarf classification is still preliminary. Three 
of the candidates are identified as quasars at z > 5.7: SDSS J141111. 29+121737.4 (z = 5.93, 
SDSS J0002+2550 hereafter), SDSS J160254.18+422822.9 (z = 6.07, SDSS J1602+4228), and 
SDSS J162331.81+311200.5 (z = 6.22, SDSS J1623+3112). The discovery spectra of the first two 
objects were obtained using DIS on the ARC 3.5m in April 2003, and the discovery spectrum of 
the last object was obtained in June 2003 at MMT using the Red Channel spectrograph. We have 
subsequently obtained longer exposures of these three quasars using the Red Channel on MMT 
and MARS on the Kitt Peak 4-meter, with total exposure times of 2 - 4 hours each. These spectra 
are shown in Figure 2. As discussed above, SDSS J000552. 34-000655. 8 (z = 5.85) was discovered 
in the deep imaging of the SDSS Southern Survey region. Figure 2 shows the discovery spectrum, 
a 40-min exposure obtained with the Clay telescope and B&C spectrograph in Oct 2002. The 
discovery spectrum of SDSS J000239.39+255034.8 (z = 5.80, §2.3) was obtained at the MMT 
using the Red Channel Spectrograph in Nov 2003. This one hour exposure is shown in Figure 2. 

The finding charts of the five new quasars are presented in Figure 1. The spectra are 
flux-calibrated to match the observed z band photometry. Table 1 presents the photometric 
properties of the new quasars, and Table 2 presents the measurements of their continuum 
properties. Following Papers I and II, the quantity ABi2so is defined as the AB magnitude of 
the continuum at rest-frame 1280A, after correcting for interstellar extinction using the map 



14 The naming convention for SDSS sources is SDSS JHHMMSS.SSiDDMMSS.S, and the positions are expressed 
in J2000.0 coordinates. The astrometry is accurate to better than 0.1" in each coordinate. 
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of Schlegel, Finkbeiner & Davis (1998). We extrapolate the continuum to rest-frame 1450A, 
assuming a continuum shape f u oc z^~ ' 5 , to calculate AB1450. None of the five quasars is detected 
in the FIRST (Becker, White & Helfand 1995) or NVSS (Condon et al. 1998) radio surveys. The 
discovery of these five new quasars brings the total of z > 5.7 quasars known to twelve, all selected 
from SDSS imaging. 



3.1. Notes on Individual Objects 

SDSS J000239.39+255034.8 (z = 5.80). This object is selected in a non-standard SDSS 
scan in the Southern Galactic Cap. It is outside the main survey region, and the spectroscopic 
follow-up of this region for z-dropout candidates is not yet complete. We therefore do not include 
it in the luminosity function calculations below. 

SDSS J0002+2550 has i = 21.47 and z = 18.99, but is undetected in 2MASS (Skrutskie 
et al. 1997), indicating that J > 16.5 and z — J < 2.5. If it were an L dwarf, it would have a 
z — J color of 2.5 or redder, so we targeted this object as a high-redshift quasar candidate. SDSS 
J1044-0125 (z = 5.74, Fan et al. 2000) was selected in a similar manner. 

At z A b ~ 19 and Mi 450 = -27.55 (A-model), SDSS J0002+2550 is very luminous. It is 
the second brightest quasar known at z > 5.7 so far and provides an excellent target for high 
resolution spectroscopic follow-up observations. The redshift of SDSS J0002+2550 is determined 
by the locations of the peak of the Lya emission line and of the 01 1300A emission line, and is 
accurate to 0.02. The quasar has a Lya+NV emission line rest-frame equivalent width of ~ 60A, 
comparable to that of most high-redshift quasars (Fan et al. 2001b). 

SDSS J000552. 34 000655. 8 (z = 5.85). This object is selected in the SDSS Southern Survey 
region. At zab = 20.54 ± 0.10, it is the faintest z > 5.7 quasar in our sample. Although the 
discovery spectrum (Figure 2) has low S/N, the strong Lya and NV emission lines and the strong 
Lyman break are clearly visible. The redshift is determined by the locations of the Lya and NV 
emission line peaks, and is accurate to 0.02. The emission lines in this quasar appear to be quite 
narrow: the Lya and NV emission are clearly separated. We estimate a FWHM of about 1500 - 
2500 km s~ l for individual lines, with large uncertainties due to the low S/N. 

SDSS J141111. 29+121737.4 (z = 5.93). The redshift is determined by the locations of the OI 
1300A and NV 1240A lines and is accurate to 0.02. The object has a moderately strong Lya+NV 
line (rest-frame equivalent width of ~ IOOA). This object is detected in two overlapping SDSS 
runs with consistent photometry (Table 1). 

SDSS J160254.18+422822.9 (z = 6.07). The redshift of this object is determined by the peaks 
of the Lya and NV 1240A lines and is accurate to 0.02. A weak OI 1300A line is also detected in 
the spectrum. The quasar spectrum shows a number of dark patches in the Lya and Ly/3 forest, 
although the S/N of the current spectrum is not sufficient to determine whether there is a short 
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Gunn-Peterson trough in the spectrum. 

Due to its proximity on the sky to the high-redshift cluster candidate CL1603, 
SDSS J 1602+4228 serendipitously lies in a 28 ks pointed observation (sequence rp800239) 
made with the ROSAT Position Sensitive Proportional Counter (PSPC). We do not find any 
significant X-ray detection of SDSS J 1602+4228 in these data. The 3<r upper limit on its 
observed-frame, Galactic absorption-corrected, 0.5-2 keV flux is 5.9 x 10~ erg cm -2 s _1 , 
adopting a power-law model with a photon index of T = 2 and the Galactic column density of 
N u = 1.3 x 10 20 cur 2 . Given the AB 1450 magnitude of SDSS J1602+4228, comparison with 
Figure 2 of Vignali et al. (2003) shows that the X-ray upper limit is consistent with X-ray 
observations of other z > 4 quasars. The slope of a nominal power law between rest-frame 2500 A 
and 2 keV is constrained to be a ox > 1.2. 

SDSS J162331. 81+311200. 5 (z = 6.22). This is the highest redshift quasar presented in this 
paper, and is the third highest redshift quasar yet known. 15 This quasar has two striking features. 
First, it has an extremely strong Lya emission line. The total equivalent width of Lya+NV is 
^ 150A in the rest-frame, without taking into account the strong absorption due to the Lya forest 
on the blue side of the Lya emission. For comparison, Fan et al. (2001b) measured the mean and 
standard deviation of the rest-frame Lya+NV equivalent width of 69.3 + 18. OA, based on a sample 
of ~ 40 quasars at z ~ 4 selected from the SDSS (see also Schneider, Schmidt & Gunn 1991 and 
§4) . The line strength of this quasar is more than a factor of two larger than the average at high 
redshift. It is one of the strongest-lined quasars at z > 4 yet known, and the strongest-lined quasar 
at z > 5. 

Second, SDSS J1623+3112 has a complete Gunn-Peterson (1965) trough. Following Becker et 
al. (2001), Fan et al. (2002) and White et al. (2003), we define the transmitted flux ratio as: 

T{zabs) = (fu bs /fu° n ) , (1 + *abs ~ 0.1) X 1216A < A < (1 + Z abs + 0.1) X 1216A, (1) 

where f£° n is the continuum level extrapolated from the red side of Lya emission. Using 
the spectrum presented in Figure 2, we find that at z = 6.05, the transmitted flux ratio 
T{ z abs = 5.95 — 6.15) = 0.004 ± 0.008, consistent with zero flux in the Gunn-Peterson trough 
region. We also detect a complete Ly/3 Gunn-Peterson trough in this quasar. This is the fourth 
quasar with a complete Gunn-Peterson trough, after SDSS J1306+0524 (z = 6.28, Becker et al. 
2002), SDSS J1148+5251 (z = 6.42, White et al. 2003) and SDSS J1048+4637 (z = 6.18, Fan et 
al. 2003). 

In all quasars at z > 6.1, complete Gunn-Peterson troughs are detected, starting from 
z abs = 5.95 + 0.1, and extending to the highest redshift not affected by the quasar proximity effect. 



15 SDSS J104845. 05+463718. 3 was reported to have a redshift of 6.23 in Fan et al. (2003). Subsequent observations 
show that it is a Broad Absorption Line (BAL) quasar (Maiolino et al. 2004, Fan et al. in preparation) and the 
original redshift determination is biased; the best redshift estimate is z = 6.18 based on new observations. 
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The detection of a complete Gunn-Peterson trough in SDSS J1623+3112 further confirms the 
rapid transition of the ionization state of the IGM at z ~ 6 (e.g. Becker et al. 2001, Djorgovski et 
al. 2001, Fan et al. 2002, White et al. 2003; see also Songaila 2004 for a different interpretation). 
The S/N of the spectrum of SDSS J1623+3112 is considerably lower than that of the other 
three quasars, therefore the optical depth limit we are able to place is not yet very strong. In a 
subsequent paper, we will present a detailed analysis of the constraints on the evolution of IGM 
properties using the five new quasars presented in this paper. 

4. Quasar Composite Spectrum at z ~ 6 

The spectral energy distributions of luminous quasars show little evolution out to high 
redshift. There is growing evidence from emission line ratio measurements that quasar broad 
emission line regions have roughly solar or even higher metallicities at z > 4 (e.g., Hamann & 
Ferland 1993, Dietrich et al. 2003a), similar to that in low redshift quasars. Dietrich et al. (2003b) 
found the Fell/Mgll ratio to have roughly the same value in a sample of z ~ 5 quasars as at lower 
redshift, suggesting that the metallicity of quasar emission line region remains high to even earlier 
epochs. 

The sample of twelve quasars at z > 5.7 from the SDSS provides the first opportunity to 
study the evolution of quasar spectral properties at z ~ 6, less than 1 Gyr after the Big Bang and 
only 700 million years from the first star formation in the Universe (Kogut et al. 2003, Spergel et 
al. 2003). Optical and infrared spectroscopy of some z ~ 6 quasars already indicates a lack of 
evolution in the spectral properties of these luminous quasars: Pentericci et al. (2002) show that 
the CIV/NV ratio in two z ~ 6 quasars are indicative of supersolar metallicity in these systems. 
Freudling et al. (2003) and Barth et al. (2003) detected strong Fell emission in the spectra of four 
z ~ 6 SDSS quasars. In addition, the optical-to- X-ray flux ratios and X-ray continuum shapes 
show at most mild evolution from low redshift (e.g. Brandt et al. 2002, Vignali et al. 2003). These 
results, if confirmed with a larger sample, suggest that the accretion disk and photoionization 
structure of quasars reached maturity very early on and are probably insensitive to the host galaxy 
environment. 

Figure 3 shows the composite of eleven of our twelve z > 5.7 quasar spectra. We omit SDSS 
J0005-0006 due to its low S/N. To produce the composite, we simply redshift all the spectra to 
zero, scale the continuum level of each quasar based on its mi45o magnitude, and average all the 
scaled spectra with equal weighting. The composite covers rest-frame wavelengths from 11 00 A 
to 1450A. Also plotted in Fi gure 3 is the low-redshift SDSS quasar composite of Vanden Berk et 
al. (2001). The effective redshift of the low-redshift composite in this redshift range is z ~ 2. 
Blueward of Lya emission, the strong IGM absorption at z ~ 6 almost completely removes the 
quasar flux. Redward of Lya emission, there is no detectable difference in the intrinsic UV spectral 
properties. The continuum shape is consistent with the power law f v cx v~ 0A measured by Vanden 
Berk et al. (2001). Clearly, spectral coverage in the near-IR is needed to put stronger constraints 
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on the continuum shape. Pentericci et al. (2003) use IR photometry of a sample of quasars at 
z = 3.5 — 6 to measure a continuum shape of f u oc i^~ 0,5 , independent of redshift. The strengths 
of the emission lines at z ~ 6, including NV 1240A, 01 1300A, CII 1335A, and SilV+OIV 1400A, 
are also comparable to those at low redshift. This composite does not go red enough to cover the 
CIV 1549A line, so we cannot test whether the Baldwin (1977) effect exists at these redshift. The 
strength of the red wing of the Lya emission line shows no evolution from that at low redshift. 
The weaker Lya emission in the blue wing is due to the strong IGM absorption. The average 
FWHM of emission lines is ~ 6000 km s _1 , also consistent with the low-redshift average. 

5. Luminosity Function of z ~ 6 Quasars 

In Papers I and II, we estimated the comoving density of quasars at z ~ 6 using a sample 
of six quasars, covering a total area of 2870 deg 2 . We repeat the calculation here, including 
the additional three quasars in the complete sample; as explained above, SDSS 0002+2550 and 
SDSS J0005-0006 were not selected as part of the flux-limited complete sample and will not be 
included in the quasar luminosity function calculation. These nine quasars form a complete sample 
satisfying the selection criteria in Eq. (1) over a total area of 4578 deg 2 . Following Papers I and 
II, we calculate the selection function of z ~ 6 quasars using a Monte-Carlo simulation of quasar 
colors, taking into account the distribution of quasar emission line and continuum properties, Lya 
absorption, the SDSS photometric errors and Galactic extinction. The selection function as a 
function of redshift z and absolute magnitude M1450 is illustrated in Figure 7 of Paper II for the 
A-model. The total spatial density of quasars at z ~ 6 is derived using the l/V a method. We 
find that at the average redshift of (z) = 6.07, p(M U50 < -26.4) = (10.5 ± 3.9) x 10~ 10 Mpc~ 3 
for the Q = 1 model, and p{M U50 < -26.7) = (6.4 ± 2.4) x 10~ 10 Mpc~ 3 for the A-model. The 
results, which are consistent with those in Paper II with smaller error bars, are plotted in Figure 4, 
together with the measurements at lower redshifts from the 2dF survey (Boyle et al. 2000, Croom 
et al. 2004) at z < 2.2, from Schmidt et al. (1995) at 2.7 < z < 4.8 and from Fan et al. (2001a) at 
3.6 < z < 5.0. The comoving density of luminous quasars at z ~ 6 is 30 times smaller than that 
at z ~ 3. 

Following Paper II, we derive the bright-end slope from the luminosity distributions of the 
sample using a maximum likelihood estimate. Assuming a single-power law luminosity function: 

*(M 1450 ) = ^*10- a4 [ Ml450+26 ]( /3+1 ), (2) 

we find that for the A-model, = (3.3± 3 ;g) x 10~ 9 Mpc" 3 . The best-fit bright-end slope is 
[3 = —3.2, with a 68% confidence range of [-2.5, -4.0] and a 95% confidence range of [-2.2, -4.2]. 
This slope is steeper than that measured in Paper II, where the best-fit value was -2.3 with a 
68% range of [-1.6, -3.1]. The reason for this 1-cr change is that all three new quasars included in 
the sample have zab ~ 20, close to the detection limit. Their inclusion, when correcting for the 
lower completeness at the faint end, drives the best-fit luminosity function to steeper slopes. For 
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example, the average selection probability of quasars with the redshift and luminosity of SDSS 
J1623+3112 is of the order 10%; the strong Lya emission boosts the z-band magnitude and makes 
the z — J color bluer, making the selection of this object easier than an average z = 6.2 quasar. 
The best-fit slope omitting SDSS J1623+3112 is -2.8. The slopes derived here and in Paper II 
differ only at the l-a level. The errors are large because of the small number of objects in the 
sample, especially at the faint end. This underlines the need for a large sample to put strong 
constraints on the quasar luminosity function at the highest redshift. 

The slope of the quasar luminosity function has important implications. The total quantity 
of ionizing photons emitted by the high-redshift quasar population is determined by the quasar 
luminosity function. Fan et al. (2001) found that z ~ 6 quasars are not likely to be the sources 
responsible for reionizing the universe, or keeping the universe ionized at high redshift, assuming 
a slope of (3 > —3.5. Interestingly, we are approaching this limit with the current analysis. 

In a flux-limited sample, the lensing probability of the brightest observed objects is boosted 
by magnification bias (e.g. Turner et al. 1984). The theoretical prediction of the fraction 
of strongly-lensed quasars at high redshift could be be of order unity for a sufficiently steep 
luminosity function (e.g. Wyithe & Loeb 2002a, Comerford, Haiman & Schaye 2002). Fan et al. 
(2003) and Richards et al. (2004) used the lack of multiply-imaged quasars among the SDSS z ~ 6 
quasar sample to constrain the shape of the quasar luminosity function to be f3 > —4.6 at the 3-<r 
level. The lensing fraction increases by a factor of ~ 5 by assuming a slope of (3 = —3.3 rather 
than P = —2.2 (Wyithe & Loeb 2002b). The quasar emissivity and lensing fraction also depends 
strongly on the faint-end quasar slope, which is currently completely unknown. 

In this paper, we present the first quasar discovered in the faint quasar survey using co-added 
catalogs from the SDSS Southern Survey region. In Fall 2003, we also started to use co-added 
images of the SDSS Southern Survey from ~ 8 SDSS observations to select candidates. This 
overcomes the limitation of the co-added catalog which requires the object to be detected in 
a single exposure. We have generated a preliminary photometric catalog using the SExtractor 
software; the integration of the co-added imaging into the SDSS photometric pipeline is currently 
underway. Using the catalogs generated by co-added imaging, we were able to recover SDSS 
J000552. 34-000655. 8. The remaining candidate identification is still in progress and will be 
reported in the future. In the next few years, the faint quasars selected from the Southern Survey 
will be combined with the bright quasars in the Main Survey region to study the evolution of 
quasar population at z ~ 6 down to much fainter luminosities to probe the evolution of faint 
quasars. 
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Table 1. Photometric Properties of Five New z > 5.7 Quasars 



object 


redshift 


i 


z 


J 


SDSS run 


J000239.39+255034.8 


5.80 ± 0.02 


21.47 ± 0.11 


18.99 ± 0.05 


> 16.5 


4152 


J000552.34-000655.8 


5.85 ± 0.02 


23.40 ± 0.34 


20.54 ± 0.10 


19.87 ± 0.10 


multiple 


J141111.29+121737.4 


5.93 ± 0.02 


23.43 ± 0.37 


19.63 ± 0.07 


18.95 ± 0.05 


3836 






22.85 ± 0.30 


19.65 ± 0.08 




3996 


J160254.18+422822.9 


6.07 ± 0.02 


22.78 ± 0.38 


19.89 ± 0.10 


18.46 ± 0.05 


3705 


J162331. 81+311200.5 


6.22 ± 0.02 


24.52 ± 0.62 


20.09 ± 0.10 


19.15 ± 0.10 


3918 



The SDSS photometry (i, z) is reported in terms of asinh magnitudes on the AB system. The 
asinh magnitude system is defined by Lupton, Gunn & Szalay (1999); it becomes a linear scale in 
flux when the absolute value of the signal-to-noise ratio is less than about 5. In this system, zero 
flux corresponds to 24.4 and 22.8, in i, and z, respectively; larger magnitudes refer to negative flux 
values. The J magnitude is on a Vega-based system. 



Table 2. Continuum Properties of new z > 6 Quasars 



object 


redshift 


^4-Bl280 




^1450 

(f2-model) 


^1450 

(A-model) 


E(B - V) 
(Galactic) 


J000239.39+255034.8 


5.80 


19.09 


19.02 


-27.40 


-27.66 


0.037 


J000552.34-000655.8 


5.83 


20.30 


20.23 


-26.21 


-26.46 


0.033 


J141111.29+121737.4 


5.93 


20.04 


19.97 


-26.49 


-26.75 


0.025 


J160254.18+422822.9 


6.07 


19.93 


19.86 


-26.63 


-26.82 


0.014 


J162331.81+311200.5 


6.22 


20.20 


20.13 


-26.40 


-26.67 


0.022 
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Fig. 1. — SDSS z-band images of the five new z > 5.7 quasars. Each side of the finding chart is 
160". The arrow at the lower left indicates the direction of North on the finding chart; East is 90° 
counterclockwise from North. 
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Fig. 2. — Spectra of the five new quasars at z > 5.7. The spectrum of SDSS J0002+2550 is a 
60-min exposure taken with the MMT-Red Channel; the spectrum of SDSS J0005-0006 is a 40-min 
exposure taken with Magellan II and B&C spectrograph; the spectrum of SDSS J1411+1217 is 
a 120-min exposure taken with KPNO-4m and MARS; the spectrum of SDSS J1602+4228 is a 
120-min exposure taken with the MMT-Red Channel; and the spectrum of SDSS J1623+3112 is a 
co- addition of two 120-min exposures using MARS and the Red Channel. The fluxes are scaled to 
reproduce the z-band magnitude as measured by the SDSS. All spectra are binned to a dispersion 
of 5A per pixel; the spectral resolutions are between 500 and 1000, depending on the spectrograph 
used. 
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Fig. 3. — The composite spectrum of eleven z ~ 6 quasars (solid line). The spectrum of SDSS 
J0005-0002 was not included because of its low S/N. The spectrum of each quasar is redshifted, 
scaled according to its 7711450 magnitude, and averaged with equal weighting. For comparison, we 
also plot the low-redshift quasar spectral composite from Vanden Berk et al. (2001). The effective 
redshift in the 1000 - 1500A range in the Vanden Berk et al. composite is about 2. The quasar 
intrinsic spectrum redward of Lya emission shows no detectable evolution up to z ~ 6, in terms 
of both the continuum shape and emission line strengths. On the blue side of Lya emission, the 
strong IGM absorption at z ~ 6 removes most of the quasar flux. 
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Fig. 4. — The evolution of the quasar comoving spatial density at M1450 < —26.7 in the A-model. 
The filled circle represents the result from this survey. The error-bar in redshift indicates the 
redshift range covered by the z-dropout survey. The dashed and dotted lines are the best-fit models 
from Fan et al. (2001b) and Schmidt et al. (1995, SSG), respectively. The solid line is the best-fit 
model from the 2dF survey (Croom et al. 2004) at z < 2.3. 



